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Q_i Abstract : 

Qh We analyse inclusive photoproduction of polarized 3 Pi quarkonium in the framework 

_C of QCD. To separate nonperturbative and perturbative parts in the density matrix of the 

produced quarkonium we use a method , which is equivalent to the diagramatic expansion 
widely used in analysing deeply inelastic scatterings. A systematic expansion in a small 
velocity v, with which a heavy quark moves inside the quarkonium in its rest frame, is 
performed for the nonperturbative parts, and they are expressed as matrix elements in 
nonrelativistic QCD. At the leading order of v there are four matrix elements representing 
nonperturbative physics. The perturbative parts are calculated at the leading order of 
coupling constants. Some numerical results, especially, numerical results for HERA are 
given. 



1. Introduction 



A quarkonium mainly consists of a heavy quark Q and its antiquark Q. The quark 
Q or the antiquark Q carries roughly the half of the momentum of the quarkonium. In 
a quarkonium rest frame the quark Q or the antiquark Q moves with a small velocity 
v. Production of a quarkonium in high energy processes can be considered as a two- 
step process, at first, a QQ pair is produced, then a transition of the QQ pair into the 
quarkonium takes place. The production of the QQ pair can be handled perturbatively in 
QCD because the mass M of Q is large, while the transition is a nonperturbative process. 
Using an expansion in the small velocity v the nonperturbative process can be described 
by some constants defined in QCD as matrix elements. 

Earlier calculations of quarkonium production have been made only within models. 
One of them is the color-singlet model. In this model the basic assumptions are that 
only a color-singlet QQ pair can be transimitted into a quarkonium and the transition is 
represented by a wavefunction of the color-singlet QQ pair. With these assumptions an 
important fact in relativistic quantum field theory is neglected, namely a quarkonium as 
a bound state consists not only of a color-singlet QQ pair, but also of many other states 
as components. Neglecting this fact leads to that there are infrared singularities in the 
perturhative results related to a P-wave quarkonium. These singularities are well known 
as the divergence in the limit of the zero-binding energy [1]. It is pointed out that a color- 
octet QQ pair combined with a soft gluon is a component of a P-wave quarkonium and this 
component is as important as a color-sinlet QQ pair [2]. A factorized formula for production 
rate of an unpolarized P-wave quarkonium is given [3], where the nonperturbative transition 
of a color-singlet and a color-octet QQ pair into the P-wave quarkonium are described at 
the leading order of v respectively by two matrix elements defined in NRQCD. With the 
factorized formula it is already shown that there is no infrared singularity in the case of 
gluon fragmentation[4,5]. 

In this work we consider inclusive photoproduction of single 3 Pi quarkonium, where 
the quarkonium is polarized. The polarization is described by a density matrix. We 
calculate the density matrix at leading order of the coupling constants and of v. To separate 
nonperturbative effects from perturbative effects we use a method, which is equivalent to 
the diagramatic expansion widely used in deeply inelastic processes[6-10]. The expansion 
provides a powerful tool to analyse and to account higher twist contributions in a deeply 
inelastic process. Although we work here at the leading order of v, one can extend with 
the method our results to higher orders. Some of our results are universal, in the sense 
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that they do not depend on the processes and the quarkonium considered here. Our 
results for the density matrix show that instead of two for unpolarized case there are four 
nonperturbative constants for polarized case, three of them are for the color-octet QQ pair. 
For unpolarized case we rediscover the factorized formula given in [3] . 

We will denote a 3 Pi quarkonium as \i- For c-quark or 6-quark it is respectively \ci 
or Xbi- The main decay of xi 1S radiative decay, where the decay products are one photon 
and a S"-wave quarkonium. This decay mode can be used to analyse the polarization 
information from the production. We will discuss how to construct spin observables for 
this purpose. Numerical results are given for photoproduction and for HERA, where the 
photon in the initial state is obtained through equivalent photon approximation from the 
initial electron at HERA. 

Our work is organized as following: In Sect. 2 we introduce the density matrix and 
separate in the density matrix the nonperturbative part from the perturbative part. In 
Sect. 3 we systematically make the small velocity expansion for the nonperturbative part 
with NRQCD and give the leading order results. Following the way explored here it is 
easy to extend the results to higher orders. Sect. 4 contains perturbative results for density 
matrices of various partonic processes at leading order of coupling constants. We discuss 
first an extension of Furry thoerem. With this the contribution from a color-singlet QQ can 
be shown to be zero at the considered order. Some numerical results are given in Sect. 5. 
We have not performed an extensive study of different observables, only some simple 
observables are studied here. In Sect. 6 is the summary of our work. For completeness we 
give in an appendix the density matrix for the radiative decay of xi ■ 



where X denotes the unobserved states, if is a quarkonium with spin 1 and A is a hadron 
in the initial state. The momenta are indicated in the brackets and the moving direction of 
the photon is opposite to the moving direction of A. We will assume that the factorization 
holds for the initial state and work only at the leading twist order, i.e., we will use the 
parton model for the initial hadron A. Therefore, we only need to study the following 
partonic process: 



2. The density Matrix and the Factorization 



We consider the following inclusive process: 



1 (q) + A(p)^H(P H ) + X 



(2.1) 



7(q) + a( qi )^H(P H )+X 



(2.2) 
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where a stands for a parton with qi = xp. In (2.2) the polarization information of the 
quarkonium is contained in a density matrix R } which can be defined as: 



R(X, X',q, qi , P H ) = Z^W? + qi ~ Ph - Px) 

x 



<H(P H ,X),X\T\ 1 (q),a(q 1 )> 
<H(P H ,X'),X\T\ 7 (q),a( qi )>* 



(2.3) 



Here T is the transition operator, A and A' denote the helicities of H . ^ means the 
average of the initial state. In (2.3) the density matrix is given in the helicity basis. Since 
helicity is not conserved, the density matrix is not diagonal in the helicity basis. Another 
useful basis is the Cartesian basis defined in the quarkonium rest frame. Throughout our 
work we define the quarkonium rest frame which is related to the quarkonium moving 
frame in (2.1) only through a Lorentz boost without rotations. Denoting the polarization 
vector as e(A) in the rest frame the density matrix Rij in the Cartesian basis is related to 
R(X, A') as 

R(X,X',q,q 1 ,p H ) = e f l {X)R lJ {q,q u p H )e J {X'). (2.4) 

Taking the Lagrangian of QED and QCD, and using Wick Thoerem one can evaluate 
the amplitude. At certain leading order of the coupling constants the amplitude takes the 
form: 

4 



<H(P H ,X),X\T\ 1 (q),a(q)> = J ^^(p^h) 



dx 4 e~^ x < H(P H ,X),X N \Q z (x)Q J (0)\0 > (2 ' 5) 
_l_ ... 

where Q{x) is the Dirac field for the heavy quark Q } the indices i,j contain the Dirac and 
color indices. The operators in the matrix element are norm ordered. The matrix element 
is a nonperturbative object, while the quantity Aij(pi,k\) is calculated by perturbative 
theory. We divided in (2.5) the unobserved state X as X = Xp + Xjy. The state Xjy is 
generated nonperturbatively and Xp is generated perturbatively with a total momentum 
as k\ . The quantity Aij (p\ , k\ ) can be obtained simply by using usual Feynman rule, and it 
can be interpreted as the amplitude for 7(?)+a(?i ) — > Q*(pi ) + Q*(<Z+<Zi — Pi ~ ki)-\-Xp(ki ). 
The real amplitude for on-shell Q and Q is uAv where the on-shell condition is taken. The 
• • • in (2.5) denotes the other possible terms by using Wick theorem. For example, there are 
terms where the matrix element contains only one operator. These terms can be regarded 
in some circumstances as contributions from parton fragmentations. The fragmentations 
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can be determined by perturbative theory [ 4,5 and references therein], they will eventually 
take after applying the procedure here for (2.5) the form as in (2.5) with the matrix element 
containing the product Qi(x)Qj(0). We will only take the contribution indicated explicitly 
in (2.5) for the density matrix in (2.3). Substituting (2.5) into (2.3) and using translational 
transformation to eliminate the summation over Xjy we obtain: 



R(\ } \' } q } q u P H ) = J2'j2 I 

X P J 



d pi o I Pi 



(2tt) 4 (2tt) 4 
Aij (Pi , h ) (To (pi ,h ho )) kl 

d 4 xd 4 x'd 4 ye- iy(q+ql - kl) - ipiX+l ^ x ' 



(2.6) 



< U\Qk{-\y)Qi{x' - \y)<f H {\')a H {\)Qi{x + \y)Qj{\y)\U > 



where a H is the create operator for H . With the form in (2.6) the contributions to the 
density matrix can be presented by the Feynman diagram in Fig.l. In Fig.l the black box 
is corresponding to the nonperturbative part: 



Tkl,ij(pi, h) = d xd x d ye 



-iy(q+ql — k 1 ) — ip 1 x + ip' 1 x' 



< n\Qk(-\y)Qi(x' - \y)<f H {\')a H {\)Qi{x + \y)Qj{\y)\U >, 



(2.7) 



while the blank part is a perturbative part and can be obtained by usual Feynman rules 
except the two pair of quark-lines connecting to the black box. These lines only indicate 
momentum flow and certain contractions of Dirac and color indices. The diagram in Fig.l 
is similar to those appearing in using diagramatic expansion for analysing deeply inelastic 
processes. 

The expression in (2.6) can be simplified by decomposing the color and Dirac indices. 
The decomposition of the color indices can be easily done with SU(3) color-symmetry. We 
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obtain the following: 

R(X, X',q, q u P H ) = R (1) (\ X',q, qi,P H ) + R (8) (\ A', q, q u P H ), 
PS 1) {X,X',q,q 1 ,P H ) = J 4^45r Trc °l° rA ^'(-P 1 ' A;i ) 

■ Tr co i or (7 A t (pi, h ^ki ■ r^(pi, h), 



Tr color T a ( 7o At( Pl , k^u ■ T^jip^h) 

^ (r>, IcA — ^-r^-v' rl 4 -iip~ i y( q+ql ~ klS >~ ipiX + i P'i x ' 



r^^(pi,A;i) = J d xd x d ye 

< ^\Qk(-\y)Qi(x' - l -y)a^(X')a H (X)Q l (x + 1^(^)10 >, 
Y%\ j {p 1 ,k 1 ) = j d 4 xd 4 x'd 4 ye- iy(q+ql - kl) - ipiX+ip 'i x ' 

< 0\Q k (-\v)T h Qi(x' - l -y)a^(X')a H (X)Q l (x + l -y)T b Q ^-y)^ > . 

(2.8) 

Here T a (a = l,---8) are the SU(3) color-matrices. The indices i,j,k and / now are 
only Dirac indices. In (2.8) there is a contribution R^ 8 \X } A', q, q\ , Pjj ) from a color-octet 
QQ pair and also a contribution R^\X, A', q, qi } Ph) from a color-singlet QQ. We will call 
them color-octet and color-singlet components respectively. The quantity (pi, k\) and 

(8) . ... . '- 

T kl ij(pi , ki) describe the nonperturbative transition from a color-singlet QQ and from a 
color-octet QQ pair into the quarkonium. In the color-singlet model it is assumed that 

ri 8 ^(Pi,^) = o. 

We will make the small velocity expansion for I^ 1 -* and I^ 8 -* before decomposing the 
Dirac indices. After the expansion the decomposition can be easily performed. To make 
the expansion we note the fact discussed in the beginning of the introduction, this fact 
leads to that the dominant space-time dependence of the matrix elements in (2.8) is given 
by 

< 0\Q(~y)Q(x' - ±y)J H (X')a H (X)Q(x + ^)Q(^)|0 > 

r - 1 n / , 1 1 x - 1 n / 1 1 m 
~ exp{ - i-P H {x -^V- ^y) + l 2 H ^ X + 2 V+ 2 V > S ' 

(2.9) 

In the quarkonium rest frame, the matrix element has only a weak dependence on x, x', 
and y, this dependence is controlled by the small scale Mv, while the dominant time- 
dependence is given by (2.9) and the correction to this is controlled by the scale Mv 2 . 
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The difference between the space- and time-dependences naturely leads to use NRQCD 
to make the expansion in v. The another advantage of using NRQCD is due the recent 
progress on lattice study, where nonperturbative properties of a quarkonium can be well 
studied with lattice NRQCD [11]. 



3. The Small Velocity Expansion 



In this section we will work in the quarkonium rest frame. To avoid confusions due to 
too many notations and bookings we still use the same booking for momenta as in Sect. 2. 
Finally we will express our results in a covariant form. In the expansion we will have 
space-time derivative c^, we will automatically substitute it with the covariant derivative 
-D^, since this is the plausible way to maintain the gauge invariance. Before going to the 
expansion we discuss briefly how the expansion of the original QCD Lagrangian in the 
power of v to obtain NRQCD. 

We write the Dirac field Q{x) as 

Q( x ) = e - mt Q+(x) + e iMt Q-(x), (3.1) 

as discussed in the last section, the four component field Qs r (x) and Q-(x) have a weak 
dependence on x. With the equation of the motion these fields take form: 

Q+(x) = P+h+(x) - l -^P+h + (x) + 0(v 2 ) 

2M (3.2) 
Q-(x) = P-h-(x) - l Jj-J- P _h_{x) + 0(v 2 ), 

here P± = (1 ± 7o)/2, hs r (x) and h-(x) are certain four component fields. Since there 
are always the combination P-h-(x) and Pj r hs r (x) appearing in (3.2), we use two two- 
component fields tp and \ to denote Pj r hs r (x) and P-h-(x). The Lagrangian for the 
two-component fields upto v 2 can be easily obtained, which is the Lagrangian of NRQCD: 

£nrqcd = < A ul) {) + ||)V> - X f (^o - ||)X (3-3) 

The physical interpretation for the two-component fields is: the field tp destructs a heavy 
quark while the field \ creates a heavy antiquark. -Lnrqcd has the exact spin-symmetry. 
However, higher order corrections to -Lnrqcd will violate this symmetry. 

3.1 The Color-Octet Part 
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At leading order of v we only need to substitute (3.1) and (3.2) into I^ 8 -*, where the 
space-time dependence of P-h-(x) and Pj r hs r (x) should be neglected: 



ri 8 ^(pMi) = -| J 4 




(3.4) 



Here the fields are at the space-time origin, the momentum P is (M, 0,0,0). At leading 
order of v the difference Mjj — 2M, which is at order v 2 , can also be neglected(we will 
discuss this later), i.e., we can take Pjj = 2P. Since only local operators are involved in 
(3.4), the Dirac indices can be now decomposed easily. We obtain: 



Here a^i = 1,2, 3) are Pauli matrices. The 8 functions will force the quark Q and Q in the 
amplitude A(pi,k\) to be on-shell. The physical meaning of the matrix elements is also 
clear, the one without Pauli matrices represents the transition of a 1 So QQ pair into the 
quarkonium, the other one with Pauli matrices represents the transition of a 3 Si QQ pair 
into the quarkonium. Now we take the case in which the quarkonium is a 3 Pi quarkonium. 
As already pointed the transition of a 3 Si QQ pair has a probability proportional to u 2 [2], 
while the transition of a 1 So QQ can only happen at higher order than v 2 . This can be 
easily seen with arguments from perturbative theory. Since the spin configuration of a 1 So 
QQ is not the same as the 3 Pi quarkonium, a gluon emission to restore the configuration 
is needed and such emission must be due to the interaction part which violates the spin- 
symmetry in NRQCD. This part is proportional to v 3 at least. Therefore the probability 
of the transition of a 1 So QQ pair into a 3 Pi quarkonium is at least at the order v 6 . Hence 
at leading order of v only a 3 Si QQ pair contributes. 

The matrix element < O|x^a; 1 T a ?/mj^(A')ajj(A)?/>^T a 07 2 x|O > can be simplified with 
rotation invariance. We write the creation operator for xi as: 



(2tt) 12 £ 4 ( 9 + qi -h- P H )8 4 (p'i ~ P)^(pi - P) 




(3.5) 



<(A) 



e z (\)aj. 



(3.6) 
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With the operator a\ the matrix element can be written: 

< 0\x^JT a ^ H (\ r )a H (\)^T a a l x\0 > = ^6 ZJ 6 xx ,a^ 

+ ( e *(A) ej (A')- e K A )^( A '))& (8) 



+ ( e *(A) ej (A') + #A) £i (A') - ISijS^cW 



(3.7) 



and 



a 



(8) 



I < 0\x*(TiT a i>a]ai^T a <Ti X \0 > 



& (8) = J_ < OlxViTXata,- - a]a t )^T a a jX \0 > (3-8) 
c (8) = 1 < IxViTXata; + aja, - ^aJaOV^V.xlO > 

The constants, i.e., the matrix elements in (3.8) are real because of the time-reversal 
invariance of QCD. Now we are in the position to define a covariant amplitude for gener- 
ating a color-octet 3 Si QQ pair: 

4 8) (A) = TrA(P, ^)T a 7 • e *(A) 7 ' ^+ Mh (3.9) 

where Mjj = 2M and Pjj = 2P. The trace is taken for Dirac index and color index. With 
this amplitude the color-octet component is: 



P^(A, A', q, qi ,Ph) = Y1 S^) 4 ^ + ?i - - P H) 

x P 



A" 

+ 6< 8 > [T^ 8) (A)T^ 8) *(A') - t£\\')T ( h 8) *(\)] (3 - 10) 
+ [T^ 8) (A)T^ 8) *(A') + T^ 8) (A')T^ 8) *(A) 

-|^Ei t h 8) ( a ")i 2 ]} 

A" 

This form is a covariant form. The summation over the color index a in (3.10) is un- 
derstood. There are three nonperturbative constants as NRQCD matrix elements in the 
color-octet component of the density matrix. No simple relation can be found between 
them. To go beyond the leading order one should simply take the space-time dependence 
of quark fields into account, except this one should also include the correction from the 
binding energy to obtain correct kinematic. At the leading order we have taken Mh = 2M. 
For the quarkonium with energy Eh, the binding energy leads to a correction proportional 
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to MfjV 2 / Efj. It is interesting to note that this effect may be included by substituting 
M with \M H in (3.1), and one obtains then instead of NRQCD an effective theory with 
Mh as a parameter. However, with such an expansion, Q and Q will be not on-shell in 
the amplitude A, it may violate gauge invariances. We do not push the issue further and 
leave it for future study. 

After summation over the polarization only the term with a^ 8 -* leads to contributions 
to the unpolarized cross section. This is corresponding to the factorization formula given 
in [3], our a*- 8 -* defined in (3.8) is related to the matrix element defined in [3] (see notation) 
as: 

a (8) = ^ <0 |Of ( 3 Si)|0>. (3.11) 
The appearance of 2Mjj is because the relativistic normalization of the hadron state used 

(8) 

here. If parity is conserved and there is no absorptive part in the amplitude T H , the term 
with M 8 -* will not contribute to the density matrix. 



3.2 The Color Singlet Part 



To expand the color-singlet part the vacuum saturation technique will help greatly. 
In the case of quarkonia this technique is a well controlled approximation [3], with it and 
(3.1) we obtain: 

< 0\Q-k(-\y)Q+i(x' - l y )a| f (A')|0 > 

< 0\a H (\)Q +t (x + I^Q.^I^IO > -{1 + 0(v 4 )} 

(3.12) 

With (3.12) we only need to analyse matrix elments containing two quark fields. Performing 
a translational transformation we have: 

< 0\a H (\)Q +t (x + l -y)Q_ J ( l -y)\Q > = < Q\a H {\)Q + i{\x)Q-j{-^x)\Q > -{1 + 0(v 2 )} 

(3.13) 

The correction in (3.13) is due to the binding energy. Here we can see that at leading order 
of v, our method is equivalent to the method used in the color-singlet model, where one 
uses wavefunction to project out different states of a color-singlet QQ pair. The analysis of 
the matrix element in (3.13) is straightforward. For H being a 3 Pi quarkonium we expand 
the matrix element in the r.h.s. of (3.13) by using (3.2) and use the parity-conservation of 
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QCD. We obtain at the leading order of v: 

- 1 1 i ( i ^ \ 

< 0\a H (X)Q+i(-x)Q-j(--x)\0 >= -— < 0|a H V' t f -- D l2 a ls J x|0 > e lll2l3 

*n°ntU X ) £ nnj 3 <> / „ M\ -^^(A) J.. 

(3.14) 

where £i 1 i 2 i s is the total antisymmetric symmbol in three dimensions with £123 = 1. 

To give our results for the color-singlet component we define a nonperturbative con- 
stant a^ 1 -* and a covariant amplitude for generating a color-singlet 3 Pi QQ pair 



= ^1 < 0|a H V> f D l2 x|0 > e Jlj2j3 | 2 , 



(3.15) 

• Tr hWfr ■ P H + M H ) ~ J^ A (P, *i)7"7 ■ Ph^ 

A{P } k l} - \ P1=P . 

With these notataions the color-singlet component of the density matrix is given by 

i?W(A, X' ,q, qi , P H ) = Y! I](2vr) 4 ^ 4 (g + ?i " h - P H ) 

x P (3.16) 

aV[T%\\){T%\\')T] 

This result is equivalent to the result in the color-singlet model(See [12]), only the non- 
perturbative effect is expressed with the NRQCD matrix element a^ 1 -* instead of a wave- 
function. In contrast to the color-octet component we have here only one nonperturbative 
constant, with the vacuum saturation it is related to the operator given in [3] through: 

a(1) = 2Mh <0 | O xi ( 3 Pi) | A0> (317) 
o 

The constant is proportional to v 2 because Di appears twice in the color-sinlet 
component is then at order v 2 . As we have seen that the color-octet component is also 
at order v 2 , hence we need to include both components in a systematical expansion of 
v. With the results in (3.16) and in (3.10), the perturbative effect is in the amplitude 

and while the nonperturbative effect is included in the four matrix elements 

defined within NRQCD. These matrix elements should be regarded as defined at a certain 
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renormalization scale, a renormalization group equation for a^ 8 -* was derived at one-loop 
level in [2] and also in [4]. 

It should be pointed out that similar method was used to analyse quarkonium de- 
cay [13], where one only takes the color-singlet QQ and the factorization was performed in 
the decay amplitude. This is reasonable since in this case the vacuum saturation is a good 
approximation. However, for the color-octet component the analysis can not be performed 
only with the amplitude, one needs to consider the density matrix i.e., the probability. 

4. The Perturbative Part 

In this section we will show that xi quarkonium can only be produced in the process 
in (2.2) at least at the order of aa 2 s . Especially, if we neglect heavy quark contents in the 
initial hadron, only the color-octet components are nonzero at order aa 2 s . Before going 
to calculate the perturbative part, we first discuss an extension of Furry theorem. This 
extension will enable us to decide which Feynman diagrams to be calculated. Consider 
a 3 Lj quarkonium to be produced via the two diagrams in Fig. 2. The only difference 
between the two diagrams is the direction of the quark line. There are n photons attaching 
to the quark line, the zth photon carries a momentum ki flowing into the quark loop. 
These photons may be off-shell or on-shell. Here, only the color-singlet component is 
involved. The black circle represents the transition of the QQ pair into the quarkonium. 
For convenience we use a wavefunction tpLM to describe the transition. Neglecting terms 
at the order higher than v, the contribution from the left diagram in the quarkonium rest 
frame can be written: 

AL = J (^( 2 -)%°)^M(q){ (pi _J )2 _ M2 • (pi _ h _\ 2)2 _ M2 

1 

' ' {pi — k\ — k 2 ■ ■ ■ k n -i) 2 — M 2 
■ Tr[7 /tl (7 • ( Pl - h) + M)7" 2 (7 • ( Pl - h - k 2 ) + M) 

•••(7-(Pi ~ki k ni ) + M)Y n (-c(Pi -PH) + M)j-e*(S z )(j- Pl + M)] } 

(4.1) 

where e*(S z ) is the polarization vector for the total spin of the QQ pair and pi = \Ph + 
To obtain contribution at the leading order of v one only needs to expand in q the terms 
in the bracket {• • •} in (4.1) and then to perform the integral. One can show that for 
L = O(S'-wave) the leading order contribution is equal to that from the right diagram 
with a multiplier factor ( — l) 1+n , while for L = l(P-wave) the factor is ( — l) 2+n . We 
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can conclude that for L = 0(1) we do not need to consider the diagrams attached with 
even(odd) number photons and for the diagrams with odd(even) number photons we only 
need to calculate the half of them. 

Now we consider possible partonic processes in (2.2) which contribute to the color- 
singlet component of the corresponding density matrix. At leading order of coupling 
constants, i.e., at eg 2 in the amplitude, there is one partonic process 7 + G — > \i + G 
via the type of the diagrams in Fig. 3, where a QQ color singlet may be produced. In 
the diagrams the two gluons must have same color, the diagrams are hence QED-like for 
the color-singlet component. From the above extension of Furry thoerem, we know that 
the sum of the contributions from the diagrams is zero at leading order of v. Therefore 
the color-sinlet component is zero at this order. In the following we will only consider 
color-octet component. 

The partonic processes below will give at order aa 2 nonzero contribution to the cor- 
responding color-octet components: 

i(q) + q(qi) ^ xi(Ph) + q(h) (4.2) 
i(q) + q(qi) ^ xi(Ph) + q(h) (4.3) 

7(q) + G(qi)^Xi(PH) + G(k 1 ) (4.4) 

For these processes we define the usual invariant quantities: 

s = {q + qi)\ i=( qi -P H ) 2 . (4.5) 

The process (4.2) is described by the Feynman diagrams in Fig. 4. There are 6 diagrams 
for the process (4.4), one of them is given in Fig. 3, the other five can be obtained through 
permutation. At first look the diagrams containing the three gluon vertex may lead to 
generating a color-octet QQ pair, but with the extension of Furry theorem, the sum of 
these diagrams are zero at leading order of v. Further, for the 6 diagrams we need only to 
calculate half of them and the corresponding amplitude is proportional to d a b c . We express 
our results in the Cartesian basis and define three tensors: 

Cuj = q l q J ~ I 

ri ii ^~ s I „' 1 2 

2 

Cszj = qWij + q'jq'u - g<Wi • (4.6) 



Mjy ( { P°h ^ P H q q° \ 
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Note q' and q^ are momenta in the quarkonium rest frame corresponding to the momenta 
q and qi. With these notations the density matrix for the process (4.2) reads: 

K = ^ hq = (2-)%? + q °-p° H - *?) - — - - — 



%J 13 v ' y ^ * n 1J 2k 96 Mfjst 

{l^a (8) [ - 2M A H + 2M 2 H (s + t) - s 2 - t 2 ] + 2M 2 H c& [C U j + 2C 2ij + C 3ij ] } , 

(4.7) 

where the nonperturbative constants a^ 8 -* and c^ 8 -* are defined in (3.8) and the subscriber 
q indicales that the initial parton is a quark. Q q is the electric charge carried by a quark q 
in the unit of the proton charge. Because of the charge conjungate symmetry the density 
matrix for the process (4.3) is the same as given in (4.7). The density matrix for the 
process (4.4) is: 

R% = <' G = (2-)%? + q° ~ P° H ~ *?) • ^ • e 2 Q 2 Q gt • f 
1 1 1 

' (s-M 2 H f ' (t-M 2 H f ' (u - M 2 H y 



{ 



-S lja ^ [s 2 (s - M 2 H f + t 2 (t - M 2 H f + u 2 (u 2 - M 2 H ) 
- 2M 2 H c^{C Uj [(s - M 2 H f + (u- M 2 H f - 2M 2 H t] 

+ C 2ij [s 2 +P] +C 3JJ 5 2 }} 

where u = Mfj — s — i. This result is similar to the result for the process 7 + G — > J /ij) + G, 
of which the total cross-section have been calculated in [14]. The differences are in the 
nonperturbative constants and in the overall color-factors. However, this similarity will be 
spoiled at higher orders. 

With the results in (4.7) and (4.8) we complete the perturbative parts in density 
matrices for the partonic processes. Before ending this section we would like to point a 
fact about xi production through gluon fusion G + G — > Xi ■> where only a color-octet 
QQ pair is produced at first. Landau- Yang theorem does not forbide this process, from 
the theorem one can only conclude that the amplitude can only be antisymmetric in the 
color-indices of the two gluons. But the amplitude is accidentally zero at the leading order 
of u[15]. At higher orders it may become nonzero. 



5. Numerical Results 



We will use the radiative decay of xi to analyse the polarization. The whole process 
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we consider is a chain-process: 



7 (g) + A(p) ^xi(Ph) + X 
Xi -> 61 + 7. 

The photon in the decay carries momentum A; in the xi res t frame. The density matrix 
/?ij(k) for the decay is given in the appendix. We used the parton model for the initial 
hadron, hence the momentum q\ in the partonic processes(4.2-4.3) is q\ = xp. Any 
observable O measured in the process (5.1) can be predicted by: 

Here f a /A( x ) 1S the parton distribution, Af is a normalization factor so that < 1 >= 1. 
An observable can be constructed with the three momenta p, Ph and k. One can also 
use instead of k the one observed in the laborotry frame. Any experimental cuts can be 
implemented into (5.2). These cuts are not only due to experiment for clean signal, but 
also are necessary in order to select the kinematical region where our perturbative results 
can be applied. For this purpose we use a simple cut. We require that the transversal 
momentum of xi to be larger than Pt- With sufficient large Pt one can be in confident 
that the employment of parton model is correct. With this cut the contribution from 
elastic/diffractive production is discarded, only inelastic production which is handled in 
this work contributes to the process (5.1). With the cut the factor Af is the total cross- 
section of unpolarized xi • We will take xi to be a charmonium Xd and the initial hadron 
A in (5.1) is a proton. For parton distributions we simply take the HO set given in [16] with 
fi = M Xcl . We take M H = M Xcl = 3.51Gev and a s (M Xcl ) = 0.25. The cut is P T = 3GeV. 
We consider here the total cross-section (T lp {xci)- The cross-section is proportional to the 
nonperturbative constant a^ 8 -* in our approximation. An estimation of this constant in 
another convention is given in [2] and in [17]. The relation to the parameter H' 8 defined in 
[2,3] is: 

« (8) = ^ ci i^(l + 0(, 2 )), (5.3) 

The parameter H' s estimated in [2,17] lies between 1.4 and 3.0 MeV. We take H' s = 2.0MeV. 
The numerical result for the total cross section in nb as a function of ^Js^ = \/ (q + p) 2 in 
GeV is given in Fig. 5. At Pt = 3GeV the dominant contribution comes from photon-gluon 
fusion, while the processes with quark or antiquark as the initial partons contribute only 
with the fraction at order of 10 -3 . 
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Now we turn to HERA. At HEAR the initial hadron is a proton. The initial photon 
can be regarded as the one emitted by the initial electron: 

e^e + 7*. (5.4) 

With the equivalent photon approximation(See detail in [18]) this photon in the remainder 
process like in the process (5.1), can be regarded as real and moving in the same direction 
as the initial electron. This is just in the same spirit as the parton model, but the photon 
distribution can be calculated perturbatively, and the distribution is different in different 
cases. We consider the case without tagging the outcoming electron, the distribution is[19] : 

ol [, , 2E(1 — z) 1 9 , 2 — z , ,1 ^,Tni. 

where E is the energy carried by the initial electron, z is the energy fraction carried by 
the photon. The total cross section a at HERA as a function of Pt is drawn in Fig. 6. The 
cross section rapidly decreases as Pt increases. At very large Pt the contribution from 
the process (4.2) and (4.3), i.e., the initial parton is a quark or antiquark is dominant, 
and this contribution can be factorized further as a real gluon is first produced and the 
quarkonium is then produced via gluon fragmentation. But at large Pt o is too small for 
HERA to observe. 

To obtain polarization information of the produced Xd a simple way may be to con- 
struct some integrated spin observables, where the momentum k of the photon is involved. 
We construct two of them at HERA 

O p = (k.pf-^ H = (k.p H f-± (5.6) 

where k = k/|k|, p is the moving direction of the proton and Ph is the moving direction of 
Xd at HERA. Special care must be taken for these observable because they are not Lorentz- 
invariant. However, it should be noted that the formula in (5.2) is Lorentz-invariant in 
the sense that < O > is invariant provided O itself is invariant. The expectation values 
of P and Oh are proportional to c^ 8 -* / a^ 8 -* and no information is available for c^K The 
numerical results for < Oh > as a function of Pt are given in Fig. 7, where we simply take 
The expectation value approaches to a constant as Pt > 8GeV. At Pt = 2GeV 
we explicitly give our numerical results in the following: 

c (8) 

a = 0.025nb, < H >= -0.012—-. (5.7) 
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The expectation value of P is much smaller than < Oh >• Hence we do not present its 
results. 

6. Summary 

In this work inclusive photoproduction of polarized 3 Pi quarkonium is treated in the 
framework of QCD. The whole information about the production is contained in density 
matrices. We have separated in these matrices the nonperturbative and perturbative parts. 
The nonperturbative parts are worked out at leading order of v, at this order there are 
four matrix elements defined in NRQCD. Three of them represent nonperturbative effects 
of the transition from a color-octet QQ into 3 Pi quarkonium, one is for the transition of a 
color-singlet QQ pair into 3 Pi quarkonium and is corresponding to the wavefunction at the 
origin in the color-singlet model. With the method employed here one can systematically 
perform the small velocity expansion and find corrections in higher order in v. Although 
we treated the case of the photoproduction, the factorized form shown here is general and 
the four matrix elements are universal. The perturbative parts are calculated at leading 
order of coupling constants. At this order they only receive contributions from color-octet 
component, i.e., a QQ pair in color-octet is first produced, then this pair is transmitted into 
a 3f quarkonium nonperturbatively. At higher order in a s the color-singlet component 
will become nonzero. 

Some numerical results are given in this work. We calculated the total cross section of 
photoproduction and also for HERA. The cross section is not very small for observation. 
Two spin observables are constructed for HERA, the numerical results for one of them is 
given. If one can measure these observables one can determine one of these four matrix 
elements, which is unknown now. Our numerical results are based on leading order results 
of the density matrices and they can be changed if one adds higher order corrections. In 
general the corrections from next leading order in a s are numerically important, also, large 
corrections from the next leading order of v for charmonium are expected. Due to lack 
of these corrections we have not performed the numerical study extensively. This can be 
done once the corrections are known. 
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Appendix 

In this appendix we consider the radiative decay of xi m its rest frame: 



3 



Si +7 



(Al) 



where the photon carries momentum k. The mass difference between the two quarkonia is 
at order of v 2 . We will assume that the polarization of the 3 S\ quarkonium and the photon 
is not observed. The density matrix pij in this case can be defined in a similar way as 
(2.3). The decay amplitude can be expanded in v. Assuming C- and P-symmetries to be 
held, at leading order the amplitude reads [20]: 



Here / is a form-factor, other notations are self-explained. With this amplitude we can 
obtained the density matrix pij at the leading order of v: 



< Si,7|T|xi >= fe^ap 



ufi *v *a U 

k e y fc3 Si e Xi . 



(A2) 




(A3) 



where k = k/|k|. The density matrix is normalized: 




(A4) 



where is the solid angle of k. 
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Figure Captions 

Fig.l: The Feynman diagram representing Eq.(2.6). The time-direction is from left to 
right. The broken line is the Cutkosky cut. The black box represents the nonperturbative 
part while the blank part is the perturbative part. 
Fig. 2 The Feynman diagrams for quarkonium production 

Fig. 3 One of the Feynman diagrams for quarkonium production through photon-gluon 
fusion. 

Fig. 4 The Feynman diagrams for process 7 + q — > QQ + q. 

Fig. 5 The total cross section in nb for Xd photoproduction as function of ^Js lp in GeV. 
Fig. 6 The total cross section in nb for HERA as a function of Pt in GeV. 
Fig. 7 The expectation value < Oh > f° r HERA as a function of Pt in GeV. 
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